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We examine a Bose-Einstein condensate of atoms that rotates in a quadratic-plus-quartic potential. 
It is shown that states of different circulation can be metastable. As a result, we demonstrate 
that the gas can exhibit hysteresis as the angular frequency of rotation of the trap is varied. The 
simplicity of the picture that emerges for small coupling strengths suggests that this system may be 
attractive for studies of phase transitions. 
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I. INTRODUCTION 

The phenomenon of hysteresis is wide-spread in physics 
and can be expected whenever a physical system under- 
goes a first-order phase transition. Consider, for exam- 
ple, a system whose free energy has two minima as a 
function of its order parameter. As some external param- 
eter is varied, these can shift with respect to each other. 
Provided only that fluctuations are not of a strength suf- 
ficient to bring the system to the absolute minimum of 
its free energy, it will remain at a local minimum even 
of this is not the absolute minimum, i.e., even if it is in 
a metastable state. The existence of metastable states 
implies that the actual state of the system depends on 
its history. 

Superfluids number among the many physical systems 
which can display hysteresis. The rigidity against weak 
perturbations of the state of a superfluid that rotates in a 
bucket or flows in a pipe can be attributed to its metasta- 
bility [1]. For this reason, these fluids are expected to 
show hysteresis and actually do so. For example, if such 
a fluid is inside a container which slowly starts to rotate, 
it remains stationary until the angular frequency reaches 
some critical value. In the reverse process, with the fluid 
in motion, the system follows a different path, and the 
fluid continues to rotate even when the container does 
not. 

In the present paper we demonstrate how these ideas 
apply to a Bose-Einstein condensate of alkali-metal 
atoms that rotates in a quadratic-plus-quartic potential. 
Our motivation for this study comes from recent experi- 
ments in such trapping potentials which observed vortex 
states [2,3]. Using a simple but still realistic model [4,5], 
we demonstrate in Sec. II the existence of metastable 
states in these systems. In Sec. Ill we examine how the 
effects of hysteresis can appear as the frequency of rota- 
tion of the trap changes. In Sec. IV we discuss our results, 
and Sec. V summarizes the basic points of our study. 



II. MODEL - STRUCTURE OF THE ENERGY 
SURFACE 

To see how hysteresis works in this problem, it is useful 
to give a brief summary of the results of Refs. [4,5]. One 
of the most important elements of the present study is to 
understand the structure of the energy of the gas in the 
rotating frame of reference, E' — E—'L-Vl, where E is the 
energy in the lab frame, L is its angular momentum, and 
rj is the frequency of rotation of the external trapping 
potential V{p). This trapping potential is assumed to 
have the form 

np) = ^M^V[i+A(^f]. (1) 

Here, p is the cylindrical polar coordinate, M is the 
atomic mass, ao = (Ti/Afw)^/^ is the oscillator length, 
and A is a small dimensionless constant. (In the experi- 
ment of Ref. [2], A ~ 10"'^.) The trapping along the axis 
of rotation, the z axis, has been neglected for simplicity, 
and the appropriate density is thus the number of atoms 
N per unit length Z,a = N/Z. 

The (repulsive) interaction between atoms is assumed 
to have the usual form 

Mnt = ^f/oE'^('-^-^^)' (2) 

where Uq = A-Kfi^a/M is the strength of the effective two- 
body interaction with a equal to the scattering length 
for atom-atom collisions. In the limit of weak interac- 
tions considered here, the typical atomic density n is 
^ N/ (■na^Z). Thus, for a typical interaction energy uUq, 
nUo/huj ^ aa and the dimensionless quantity which plays 
the role of a coupling constant in this two-dimensional 
problem is era, which is assumed to be smaller than unity. 

In the limit era ^ 1 and A 1 the appropriate basis 
states with angular momentum mh are the eigenstates of 
the harmonic potential with no radial nodes, 
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where (j) is the angle in cyUndrical polar coordinates. 
Quite generally, the order parameter, of the conden- 
sate can then be expanded in this basis as 



(4) 



It was shown in Ref. [4] that, for sufficiently small aa, the 
energy of the gas in the rotating frame, E' , is minimized 
when only one component m = toq in the expansion of 
Eq. (4) contributes to the order parameter, = . In 
addition, the critical frequencies which denote the lower 
limit on the absolute stability of the states are given 
by 
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The solid (straight) line in Fig. 1 shows this boundary 
between the states with mo = and mo = 1. On the 
left of this line the energy E' has an absolute minimum 
for = $o; on the right the absolute minimum occurs 
for '5 = $1. This is shown schematically in Fig. 2, where 
the minimum on the left in the graphs corresponds to the 
state $0, while the one on the right corresponds to the 
state 3>i. 

As discussed in Refs. [4,5], there are metastable states 
for certain values of aa and fi/w in addition to these 
absolute minima. Consider, for example, the case mo = 
of a non-rotating condensate. (The same procedure 
can be applied to any other mo.) While the state $o 
represents the absolute minimum of E' in regions I and 
II, it is metastable in region III and unstable in region 
IV with respect to the state 



* = C_i$_i + Co$0 + Ci$i 



(6) 



for aa < 0.06. For aa > 0.06 (but aa still sufficiently 
small) the instability is towards the state 



* = C_2$-2 + Co$0 + C2$2- 



(7) 



These genuine instabilities are given by the right dashed 
line in Fig. 1. 

To be specific, the energy £" of the system per particle 
in the rotating frame in the state of Eq. (6) is 
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Making use of the normalization condition, E' can be ex- 
pressed in terms of cq and ci , and is a two-dimensional en- 
ergy surface. As one crosses the line Cl/co = l + 2X — aa/2 



(i.e., the straight solid line in Fig. 1 given by Eq. (5) for 
mo = 0), the absolute minimum moves discontinuously 
from |cop = 1 to |cip — 1. Furthermore, even above 
this line, there is a region in the O/w - aa plane (i.e., 
region III in Fig. 1), for which E' has a local minimum at 
|co|^ = 1. The state $o is thus an absolute minimum of 
E' in regions I and II, and it describes a local minimum 
in E' in region III, i.e., it is metastable. In region IV, 
on the other hand, the local minimum disappears. These 
observations are shown schematically in Fig. 2, where the 
minimum on the left corresponds to $o- 

Similarly, for ^' = $i, the left dashed line in Fig. 1 
gives the boundary for the (local) stability of the state 
$1 with respect to the state 



* = Co$0 + Ci$i + C2$2- 



(9) 



In regions III and IV, <I>i represents the absolute mini- 
mum in E' , in region II it is metastable, and in region I 
it no longer represents a local minimum. These facts are 
also shown schematically in Fig. 2, where the minimum 
on the right corresponds to 

The general picture of stability/metastability that 
emerges for all mo > 2 is shown as the dashed line in 
Fig. 3 for mo = 2. The state $2 provides the abso- 
lute minimum of E' between the straight lines, while in 
the remaining part, $2 is metastable. Along the dashed 
curve, and as one moves clockwise, the instability is 
against a linear combination of states involving the fol- 
lowing m: (1, 2, 3), (0, 2, 4), (-2, 2, 6), (-1, 2, 5), (0, 2, 4), 
and (1,2,3). 



III. HYSTERESIS 

With these observations, it is easy to demonstrate the 
effect of hysteresis. Assume that the condensate is ini- 
tially at rest and that CI increases slowly from zero with 
aa kept fixed (and less than 0.06 for simplicity.) So long 
as CI/lu is in the regions I and II, the cloud will not ro- 
tate since E' has its absolute minimum in the state with 
mo = 0. Provided that fluctuations are sufficiently weak, 
the system will not rotate even in region III since <I>o is 
a metastable state. The gas will begin to rotate only 
when Q/u) is in region IV, where $0 is unstable against 
the state of Eq. (6). The system will then move to the 
pure state $1 (i.e., the absolute minimum of the system) 
on the fastest available time scale. In this process, it is 
important to note that the gas is unstable with respect 
to the state of Eq. (6) in region IV. Since $1 is one of 
the components of there is a monotonically-decreasing 
path in the energy from the (unstable) pure state $0 to 
the pure state <I>i, which is the absolute minimum of E'. 

If one now considers the inverse process in which f2 
decreases slowly from region IV with the gas initially in 
the state $1, the system will remain in this state in re- 
gions IV and III where it is the absolute minimum and 
in region II where $1 is metastable. In the absence of 
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FIG. 1. The phase diagram of a rotating Bose-Einstein con- 
densate confined in a quadratic-plus-quartic potential in the 
Q,/(jj - aa plane. The straight solid line gives the phase bound- 
ary between the phases where the states $o (regions I, II) and 
$1 (regions III, IV) provide the absolute minima of the en- 
ergy in the rotating frame. $o is mctastablc in region III and 
unstable in region IV. Similarly $i is mctastable in region II 
and unstable in region I. 

significant fluctuations, a transition will occur only when 
fl/w reaches region I. In this region, $i is unstable to 
the state of Eq. (9), which includes a component of $o- 
The system will move to the absolute minimum described 
by the (pure) state <I>o, and rotation will cease. The gas 
thus displays hysteresis in regions II and III. The order 
parameter is given by $o as increases; by $i as it 
decreases. 

A similar process evidently exists for larger values of 
f2, and we show the first few steps of the angular mo- 
mentum m of the gas as a function of O/w in Fig. 4 for 
A = 0.05 and aa = 0.05. The picture that emerges for 
small values of aa is clear, since the leading instability of 
a given state (for increasing/decreasing O) always 
involves a linear superposition of the form 

with n = 1. (For aa = the equation that yields 
each solid curve coincides with those leading to the cor- 
responding dashed curves [4,5]. Thus, the two dashed 
curves always meet their solid partner at era = 0.) As 
soon as becomes unstable, the order parameter will 
adjust rapidly to the new minimum provided by the state 

We have found that the largest value of era, (cra)c, for 
which the instability of a given ^rno is to a state of the 
form of Eq. (10) with n = I increases roughly linearly 
with m [{aa)c = 0.06 for mo — 0]. In this case E' can 
have at most two local minima. The transition between 
different states as fl increases/decreases is then clear and 
consists of steps in the angular momentum m of the gas 
equal to ±1. For fixed aa, the pattern shown in Fig. 4 
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FIG. 2. Schematic diagram of the energy of the system in 
the rotating frame as a function of the order parameter for 
a fixed aa and for values of fl/ui in regions corresponding to 
those in Fig. 1. The minimum on the left corresponds to the 
state 3>o; that on the right to the state 3>i. 

will always hold if Q/u is sufficiently large since {aa)c 
increases with m and aa will eventually become smaller 
than {aa)c- 

For values of aa > (aa)c the situation becomes more 
rich, since there are instabilities to many different com- 
binations of states, depending on the specific value of aa 
[like, for example, the one in Eq. (7)] of the general form 
given in Eq. (10), with n being a integer. In this case 
there is a wide variety of possibilities where the system 
will evolve as n changes. The (multi-dimensional) en- 
ergy E' develops many local minima and small changes 
in aa may drive the gas to completely different states, 
as it will follow the path with the steepest decrease in 
its energy. The corresponding graph of Fig. 4 is then in 
principle more complicated, as the steps in \m\ may be 
different than unity and can be non-integer. When the 
absolute minimum energy is given by Eq. (10) with all 
coefficients non-zero, the order parameter has an n-fold 
discrete rotational symmetry [5]. This observation is of 
some importance. The existence of metastable states is a 
necessary condition for hysteresis. As repeatedly noted, 
the ability to observe hysteresis also requires the absence 
of fluctuations sufficient to cause phase transitions from 
mctastablc states. The fluctuations required to change 
discrete rotational symmetry become relatively improb- 
able and the experimental indications of hysteresis are 
likely to be clearer for larger values of aa where one or 
both states are given by Eq. (10) [5]. 

IV. DISCUSSION 

Reference [6] has shown that similar effects of hystere- 
sis also appear in harmonically-trapped Bose-Einstein 
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FIG. 3. Regions of absolute stability/metastability of the 
state 3>2. Between the sohd hnes $2 provides the absolute 
minimum of E' , while in the remaining region enclosed by 
the dashed curve it is metastable. 

condensates. (See also [7].) This is no surprise since hys- 
teresis is a general feature of superfluids. These results 
are also consistent with the corresponding experimental 
observations of Rcf. [8] . 

One important aspect of our study, which also applies 
in the case of harmonic trapping, is that the critical fre- 
quencies for the creation of vortex states observed in such 
experiments depend on how the experiment is performed. 
For example, if one first cools the gas in the condensed 
phase and then rotates the trap, the picture that emerges 
is that of Fig. 4. If, on the other hand, the trap is first set 
in motion and the gas then cooled, it will end up at the 
absolute minima of E' as given by Eq. (5). (See Fig. 1 in 
Ref. [4].) 

The simple structiirc of E' for small values of aa de- 
scribed here is an interesting feature, and this system 
could prove to be an ideal laboratory for the investiga- 
tion of the effects of quantum and thermal fluctuations 
on phase transitions. It is important to mention that the 
boundaries we have calculated here are exact in the limit 
of small aa. As in other problems that involve phase 
transitions between metastable states, it is assumed that 
fluctuations are suppressed because of the energy barrier 
that separates the two states. In the presence of such an 
energy barrier, the characteristic time scale for the phase 
transition is expected to be long particiflarly when the 
two local minima have distinct discrete symmetry. When 
the barrier disappears and the metastable state becomes 
unstable, the phase transition proceeds rapidly. Here, 
we have suggested an essentially exact description of the 
phase boundaries and the corrersponding energy surface. 

The effects considered here could be useful for (i) prob- 
ing persistent currents in superfluids, (ii) providing ex- 
perimental information regarding the energy surface of 
the system, and (iii) investigating the effect of fluctua- 
tions on the phase transitions (e.g., by varying either the 
number of atoms or the temperature). 
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FIG. 4. Demonstration of hysteresis. The graph shows 
the angular momentum, rn, of the system as function of 
the frequency of the rotation of the trap, ft/o;, as fi/o; in- 
creases (solid lines) and as it decreases (dashed lines). The 
dashed linos have been displaced slightly. Here, A = 0.05 and 
aa = 0.05. 



V. SUMMARY 

In summary, we have investigated the energy of a 
Bose-Einstein condensate that rotates in a quadratic- 
plus-quartic trapping potential. The energy surface is 
elementary in the limit of small coupling with a pattern 
of absolutely stable and metastable states which can give 
rise to potentially observable persistent currents and hys- 
teresis phenomena. Finally, because of its simplicity, the 
present system can offer some insight into the physics of 
phase transitions that are driven by fluctuations. 
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